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Average  1/10  highest  vertical  acceleration  at  90%  LQA 
forward  of  transom 

Average  1/10  highest  vertical  acceleration  at  center  of 
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Maximum  breadth  over  chines 

Propeller  diameter 
Propeller  expanded  area  ratio 
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NOTATION 
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Propeller  torque  load  coefficient 
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ABSTRACT 


A  procedure  for  the  prediction  of  powering 
requirements  and  vertical  accelerations  at  the 
preliminary  design  stage  is  presented  for 
planing  hulls  with  propellers  on  Inclined  shafts. 

Envelopes  of  operating  speed  versus  wave  height 
are  developed  based  on  (1)  maximum  speed  in 
seaway  due  to  limits  of  the  prime  movers  and/or 
propellers  and  (2)  human  endurance  llmitB  due  to 
vertical  accelerations.  Complete  documentation 
of  the  computer  program  is  provided  in  the 
appendix.  Sample  computations  and  plots  are 
presented  for  a  typical  planing  craft. 

ADMINISTRATIVE  INFORMATION 

This  project  was  authorized  and  partially  funded  by  the  Naval  Sea 
Systems  Command  Detachment  Norfolk  (NAVSEADET  Norfolk)  Work  Request 
N64281  80  WR  0  0062  under  Work  Unit  1-1524-712.  Development  of  the 
program  was  also  partially  funded  by  the  David  W.  Taylor  Naval  Ship  R&D 
Center  (DTNSRDC)  Ship  Performance  and  Hydromechanics  Program  under 
Work  Unit  1-1500-104. 

INTRODUCTION 

NAVSEADET  Norfolk  requested  DTNSRDC  to  develop  a  computer  routine 
to  predict  the  operational  limits  of  planing  craft  in  a  seaway  based  on 
state-of-the-art  technology. 

The  computer  program  predicts  the  resistance,  thrust  requirements, 

and  vertical  acceleration  of  a  planing  hull  for  a  matrix  of  speeds  and 

significant  wave  heights.  It  also  estimates  the  maximum  thrust,  as  a 

function  of  speed,  which  can  be  developed  with  pre-selected  prime  movers, 

reduction  gears,  and  propellers  on  Inclined  shafts.  Speed,  wave-height 

envelopes  are  then  established  based  on  the  power  limits  of  the  propulsion 

system  and  the  endurance  limits  of  the  crew  due  to  accelerations.  This 

program  labeled  PHPRLM  is  completely  documented  in  Appendix  A,  with 

sample  input  and  output  in  Appendix  B.  A  similar  program  for  a  waterjet 

propulsion  system  is  documented  in  Reference  1.  Both  programs  augment 
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the  planing  hull  feasibility  model  program  PHFMOPT  and  should  eventually 
be  incorporated  into  the  master  program. 

^■References  are  listed  on  page  10. 


PROCEDURE 


THRUST  REQUIREMENTS 

The  calm-water,  bare-hull  resistance  R^  is  generally  derived  from 
the  synthesized  Series  62-65  resistance  curves  presented  in  Figure  9  of 
Reference  2,  However,  if  more  precise  data  are  available,  e.g.,  model 
experiments  of  the  exact  hull  or  a  similar  form,  this  resistance  data 
can  be  input  directly  for  the  matrix  of  speedB  considered.  Resistance 
of  the  appendaged  hull  R  is  approximated  using  appendage  drag  factors  IV 
either  from  Reference  3  or  direct  input  if  other  data  is  available; 

Rfl  ■  R^/r^*  Added  resistance  in  rough  water  R#w  is  predicted  from  an 
empirical  equation  recently  developed  by  a  regression  analysis  of  planing 
hull  rough-water  experimental  data, 

K.„M  -1.3  («1/3^PX)0-S  *„7  a,/vl/3r2-5 

where  F^y  ■  V/(g  -  -  - 

The  total  resistance  for  the  hull  then  is  R  ■  R  +  R  .  Thrust  deduction 

t  a  aw 

factors  1-t  from  either  Reference  3  or  direct  input  are  used  to  calculate 
the  thrust  requirement  T  ■  Rt/(l-t). 

PROPELLER  CHARACTERISTICS 

Propeller  open-water  characteristics  are  derived  as  a  function  of 
pitch  ratio  P/D,  expanded  area  ratio  EAR,  and  number  of  blades  Z  with 
coefficients  developed  from  the  Wageningen  B-Screw  Series.^  Thrust  and 
torque  coefficients  K^,  and  Kq  for  flat  face,  segmental  section  propellers 
such  as  the  Gawn-Burrill  Series^  tend  to  be  slightly  higher  than  the 
B-Screw  airfoil  section  propellers.  This  difference  can  essentially 
be  taken  into  account  by  varying  the  EAR  and/or  P/D  used  in  the  open- 
water  equations.  If  actual  open-water  data  are  available  for  the 
selected  propeller,  these  data  can  be  input  to  the  program. 

The  propeller  characteristics  in  a  cavitating  environment  are 
derived  from  the  maximum  thrust  and  torque  load  coefficients  xc  Qc 
developed  as  a  function  of  cavitation  number  at  0,7  radius  ^  for 
several  propeller  series  in  Reference  7.  Options  are  available  in  the 
program  for  tabulating  aud/or  plotting  (see  Figure  1)  the  propeller 


and  Kq  as  a  function  of  advance  coefficient  J  for  the  open-water  condition 
as  well  as  the  transition  and  fully  cavitating  regions  at  several  cavitation 
numbers  o  .  These  "transition"  curves  do  not  exactly  match  the  shape  of 
experimental  cavitation  data  since  they  represent  only  80  percent  of  the 
maximum  thrust  and  torque  lines  developed  in  Reference  7  and  are  not  faired 
into  the  point  of  or  breakdown.  The  80  percent  criteria  Is  based  on 
full-scale  trial  data  which  indicates  actual  thrust  and  torque  in  the 
transition  region  to  be  less  than  the  raaximums  derived  from  the  propeller 
series  data  —  see  Figures  5  and  6  of  Reference  7.  For  Gawn-Burrlll  type 
propellers,  the  10  percent  back  cavitation  line  shown  in  Figure  23  of 
Reference  6  is  used  as  a  design  criteria  for  adequate  blade  area.  The 
tabulated  output  is  marked  with  the  letter  C  to  represent  thrust  which 
exceed  the  10  percent  cavitation  criteria  or  with  a  *  when  80  percent 
of  maximum  thrust  is  attained.  The  thrust  at  which  10  percent  back  cavitation 
occurs  is  represented  by  a  dash  line  on  the  thrust-speed  curves  shown  in 
Figure  2. 

POWER  REQUIREMENTS 

After  the  thrust  requirements  are  computed  for  the  matrix  of  speeds 
and  wave  heights  desired,  J,  K^,  and  then  are  interpolated  from  the 
propeller  characteristic  curves  as  functions  of  thrust  loading  1^/J^  and 
0  ,  Corresponding  propeller  rpm  N,  torque  Q,  delivered  horsepower  F^, 
propulsive  efficiency  r^,  propeller  efficiency  r^,  etc.  are  then  calculated. 
Thrust  wake  factor  1-w  from  either  Reference  3  or  input  are  used,  with  the 
relative  rotative  efficiency  ri^  assumed  to  be  one,  i.e.,  torque  wake  equal 
to  thrust  wake  (torque  required  in  open  water  is  equal  to  that  required 
behind  the  ship) . 

ENGINE  TORQUE-RPM  LIMITS 

Engine  characteristics  are  input  as  an  array  of  engine  rpm  versus 

brake  power  Pfi  values  covering  the  operational  range  of  the  engines. 

Gear  loss  and  shaft  loss  constants,  generally  about  2  percent  each,  are 

used  to  estimate  the  corresponding  shaft  power  Pg  and  the  power  delivered 

at  the  propeller  PQ.  For  a  given  gear  ratio  m  the  corresponding  propeller 
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rpm  and  torque  limits  are  then  established.  If  the  gear  ratio  is  not  input, 
the  program  will  compute  an  optimum  for  maximum  speed  in  a  specified  sea  state. 

m 

®opt 

The  available  thrust  at  the  torque  limits  of  the  prime  movers  is  derived 
after  interpolation  of  the  propeller  data  to  obtain  J,  K^,  then  K^,  as 
functions  of  torque  loading  and  a  .  Similarly  the  thrust  at  maximum 

rpm  Is  obtained  after  interpolation  for  as  functions  of  J  and  a  . 

A  sample  illustration  of  thrust  requirements  and  thrust  limits  due  to 
maximum  rpm  and  torque  of  the  prime  movers  is  shown  in  Figure  2  for  a 
typical  planing  hull  with  diesel  engines.  Diesel  engines  generally 
operate  with  nearly  constant  torque  whereas  gas  turbines  have  nearly 
constant  rpm.  The  upper  bound  of  the  rpm  limit  curve  corresponds 
to  the  "transition"  region  of  the  propeller  characteristic  curves  and 
represents  the  maximum  thrust  which  can  be  developed  by  the  propellers 
unless  operating  in  the  supercavitating  regime,  at  low  J's. 

MAXIMUM  SPEED 

The  maximum  speed  attainable  in  a  given  sea  state  Is  represented  by 
the  intersection  of  the  thrust  requirement  curve  and  the  thruBt  limit 
curve  due  to  either  engine  torque  or  rpm  restriction,  whichever  occurs 
first.  With  an  optimum  gear  ratio,  maximum  speed  for  a  specified  sea  state 
is  attained  at  maximum  power  of  the  prime  movers,  i.e.,  intersection  of 
1  torque  limit  and  rpm  limit  curves.  The  maximum  speed  is  derived  for  each 

wave  height.  Figure  3  shows  a  sample  graph  of  wave  height  ve  maximum  spaed. 

I  The  computation*  and  plots  can  be  made  for  1  to  4  different  gear  ratios  for 

each  caee,  one  of  which  cen  be  the  optimum  m  derived  in  the  program. 

HABITABILITY  LIMITS 

u 

Average  1/10  highest  vertical  accelerations  at  the  center  of  gravity 
aCG  and  at  the  bow  *B0W  (90X  of  the  overall  length  forward  of  the 
transom)  are  estimated  for  the  matrix  of  speeds  and  wave  heights  with 

g 

empirical  equations  recently  derived  from  experimental  data. 


_ engine  rpm  at  maximum  power  of  engines 

propeller  rpm  required  in  specified  sea,  at  max,  engine  power 


4 


0.75 


aCG  “  7,0  ^H1/3^BPX^  ^  +  T/?2*  <LP7lW  FnV 

aBOW  ’  10,5  (Hl/3/BPXJ  (1  +  XtT>  (LP/BPX5  FnV 

The  accelerations  at  any  location  between  the  CG  and  the  bow,  such  as  the 
helmsman's  station,  can  then  be  approximated  by  linear  interpolation. 

These  data  are  interpolated  to  obtain  speed,  wave-height  envelopes  for 
average  1/10  highest  vertical  accelerations  of  1.0  g,  representing  the 
endurance  limit  of  the  crew  for  4  to  8  hours,  and  1,5  g,  representing 
the  endurance  limit  for  1  to  2  hours.  Sample  speed,  wave-height  envelopes 
are  shown  in  Figure  3. 

PROPELLER  SELECTION 

The  program  can  be  run  with  any  number  of  propeller  sets  to  aide 

in  the  selection  of  optimum  parameters.  P/D,  EAR,  and  Z  must  be  input 

for  each  set.  Propeller  diameter  D  may  either  be  input  or  selected  by 

the  following  method.  Design  speed  and  wave  height  are  input  and  the 

corresponding  thrust  requirement  is  computed.  A  design  power,  not 

exceeding  the  maximum  power  of  the  prime  movers,  is  also  input.  A  minimum 

diameter  D  is  computed  from  the  thrust  loading  K_/J  corresponding  to 

the  10  percent  back  cavitation  criteria  from  Reference  6.  A  maximum 

diameter  D  is  computed  from  K_/J  at  peak  open-water  efficiency.  The 
max  • 

power  and  rpm  requirements  at  design  speed  are  also  computed  for  several 
D' s  in  5-inch  increments  from  Dm^n  to  Dmx.  ®min  8e^ecte<*  as  t*le  optimum 
diameter  if  the  power  required  with  Dmin  does  not  exceed  the  design  power. 
Optimum  diameter  is  interpolated  from  the  5-inch  increment  array  of  D'r 
if  design  power  falls  between  the  power  requirements  for  D^^  and 
When  design  power  exceeds  the  requirements  for  both  Dmin  and  °max’  Dmax 
is  selected  as  the  diameter  except  for  the  following  case.  When  Dmln  is 
actually  greater  than  Dmax,  indicating  more  than  10  percent  cavitation 
at  peak  efficiency,  Dmin  is  the  diameter  selected.  The  optimum  diametur 
is  always  rounded  up  to  the  next  full  inch. 
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No  selection  process  for  P/D  or  EAR  is  built  into  this  program 
since  many  different  factors  may  influence  a  particular  design.  However, 
numerous  jets  of  parameters  can  be  run  for  each  case  at  minimal  cost 
for  use  in  the  development  of  design  charts.  In  the  diameter  selection 
routine,  diameters  corresponding  to  even  increments  of  100  rpm  are 
interpolated  from  the  array  of  5-inch  increment  diameters.  These  values 
can  be  used  for  plotting  contours  of  rpm  on  the  propeller  design  charts 
to  aide  in  the  selection  ot  a  propeller  to  match  specific  engine  and/or 
gear  ratio  requirements. 

COMMENTS 

This  program,  together  with  the  planing  hull  feasibility  model 
program  PHFMOPT ,  is  quite  useful  for  making  timely  preliminary  design 
studies  for  planing  craft.  It  can  also  be  used  for  displacement  ships 
since  the  program  has  options  for  input  of  the  resistance  and  propulsion 
coefficients,  but  thia  program  does  not  take  into  account  any  difference 
in  thrust  and  torque  wake. 

Some  means  of  predicting  roll  at  the  preliminary  design  stage  is 
also  desirable  for  enhancement  of  this  program.  A  low'-speed  roll 
criteria  neede  to  be  established  for  completely  defining  the  speed, 
wave-height  envelope. 
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APPENDIX  A 

DOCUMENTATION  OF  SUBPROGRAMS  FOR  PHPRLM 


If 


NAME:  PROGRAM  PHPRLM 

PURPOSE:  (1)  Predict  the  resistance,  thrust  requirements,  and 

vertioal  accelerations  of  a  planing  hull  for  a  matrix  of 
speeds  and  significant  wave  heights. 

(2)  Determine  maximum  speed  obtainable  at  eaoh  wave 
height  with  pre-seleoted  prime  mover(s)  and  propeller(s) . 

(3)  Determine  wave  heights  versus  speed  corresponding 


to  human  acceleration  endurance  limits 

of  1.0 

g  and  1 . 5  g 

SUBPROGRAMS: 

PRCOEF ,  PHRES,  SAVIT,  PRCHAR,  PRINTP,  HPLOT, 
TP LOT,  MINP,  YINTE,  YINTX,  CALCOMP  Routines 

PLOT  &  PLOTS 

INPUT: 

Via  Punohed  Cards 

Card 

Columns 

N CASES 

Number  of  cases— repeat  following 
set  of  cards  for  eaoh  case 

1 

1-4 

TITLE 

Identification  for  hull  design 

2 

1-80 

XLOA 

Overall  length  of  hull  (Lqa)  in  ft 

3 

1-8 

PL 

Projected  ohine  length  of  hull  (Lp) 
in  ft 

3 

9-16 

BPX 

Maximum  beam  over  ohines  (BDV)  in  ft 

3 

17-24 

HT 

Draft  at  transom  (T^)  in  ft 

3 

25-32 

DLBS 

Displacement  at  rest  (  A  )  in  lb 

3 

33-40 

BETAM 

Deadrise  angle  at  midships  (8m)  in  deg 

3 

41-48 

XLCG 

Distanoe  of  oenter  of  gravity  forward 
of  transom  (ffS)  in  ft 

3 

49-56 

VCG 

Distanoe  of  oenter  of  gravity  above 
baseline  (KG)  in  ft 

3 

57-64 

XACC 

Distance  forward  of  transom  in  ft 
at  whioh  accelerations  are  to  be 
computed,  in  addition  to  CG  and 
bow  locations 

3 

65-72 

FT  RIM 

Fixed  trim  angle  (  t  )  in  deg 

If  not  input,  program  will  use 

Savitsky  prediction  of  trim  at 
each  speed 

3 

73-80 

PEMAX 

Maximum  brake  horsepower  of  eaoh 
prime  mover  ( P^^) 

1-8 
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PROGRAM 


PHPRLM 


Card  Column! 


REMAX  Maximum  revolutions  per  minute  (rpm)  4 

of  prime  mover  ^emty) ,  corresponding 
to  ^emax 

GL  Geer  loss  ratio  (K-)  •  brake  4 

horsepower/ shaft  horsepower 

SL  Shaft  loss  ratio  (K«)  ■  shaft  4 

horsepower/hp2  developed  at  prop 

2  4 

RHO  Water  density  (p)  in  lb  x  sec  /ft  5 

VIS 5  Kinematic  viscosity  of  water  (v)  5 

in  ft2/ sec  x  10* 

GA  Acceleration  of  gravity  (g)  in  ft/sec2  5 

DCF  Correlation  allowance  resistance  5 

predictions,  generally  sero 

SDF  Standard  deviation  factor  for  5 

resistance  data.  Use  sero  for  mean 
Series  62-65  curves 

NV  Humber  of  speeds  —  maximum  of  20  6 

NWH  Number  of  wave  heights,  including  sero  6 

—  maximum  of  5 

NPR  Number  of  prime  movers  -  6 

number  of  propellers  (n^) 

NE  Number  of  points  input  from  engine  6 

characteristics  curve  —  maximum  of  10 

IPROP  Control  for  type  of  propellers  6 

1  for  Gawn-Burrill  type 

(flat  face,  segmental  sections) 

2  for  Nawton-Rader  type 

3  for  Wageningen  B  -  Screw  type 

(airfoil  sections) 

IPM  Control  for  type  of  prime  movers  6 

1  for  diesel  engines 

2  for  ges  turbines 


9-16 

17-24 

25-32 

1-8 

9-16 

17-24 

25-32 

33-40 

1-4 

5-8 

9-12 

13-16 

17-20 


21-24 
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PROGRAM  PHPRLM 


IPLOT  Control  for  graphioal  output,  i.e., 

CALCOMP  plots 
0  for  no  plots 

1  for  plots  of  thrust  va  speed  and 

speed-wave  height  envelopes 

2  for  plots  of  propeller  oharaoteristios 

in  addition  to  plots  above 

IPC  Control  for  propulsion  coefficients 

1  if  thrust  deduction  faotor  (1-t), 

thrust  wake  faotor  (1-w),  and 
appendage  drag  faotor  (na)  are 
estimated  by  Subroutine  PRCOEF 

2  if  1-t,  1-w,  and  na  are  input 

IRES  Control  for  bare  hull  resistance  data 

1  if  resistance  is  estimated  by 

Subroutine  PHRES  based  on 
Series  62-65  data 

2  if  resistance  is  input 

I0W  Control  for  propeller  open-water 

oharaoteristios 

1  if  oomputod  by  Subroutine  OWKTQ 

based  on  Wageningen  B-sorew 
ooeffioients 

2  if  points  from  open-water  ourves 

are  input 


to  omit  engine  characteristics  and 
propeller  characteristics  (pages  2,4,5) 


Array  of  ship  speeds  (Vg)  in  knots, 
in  ascending  order,  Maximum  of  20. 

Do  not  inolude  zero  speed. 

10  speeds  per  oard.  Use  2  cards  if  NV>10 


Columns 

25-28 


29-32 


33-36 


37-40 


NJ 

Number  of  points  input  from  each  open- 
water  aurve,  if  IOW  =  2.  Maximum  of  60. 

6 

41-44 

IPRT 

Control  for  amount  of  printed  output 

0  for  complete  output  (pages  1-10) 

6 

45-48 

PROGRAM  PHPRLM 


Card 

Columns 

H13 

Array  of  significant  wave  heights  (H1/3) 
in  ft,  in  ascending  order.  Maximum  of  5. 
First  wave  height  must  be  zero. 

8 

1-8 

9T16 

PE 

Array  of  operating  horsepower  values  for 
each  prime  mover  (Pe),  in  ascending 
order.  Maximum  of  10. 

9 

1-8 

9-16 

RE 

Array  of  rpm  values  for  prime  mover 
(Ne),  corresponding  to  Pa  values  on 
oard  9 

10 

1-8 

9-16 

TDF 

Array  of  1-t  values  corresponding  to 
speeds  on  Card(s)  7*  Use  2  cards  if 

NV  >  10 .  Omit  Card(s)  11  if  IPC  =  1 

11 

1-0 

9716 

1 

TWF 

Array  of  1-w  values  at  eaoh  speed 

Omit  Card(s)  12  if  IPC  a  1 

12 

1-8 

9-16 

ADF 

Array  of  values  at  eaoh  speed 

Omit  Card(s)  13  if  IPC  =  1 

13 

1-8 

9-16 

RBH 

Array  of  bare  hull  resistance  (Rj,) 
in  lb  at  eaoh  speed.  Appendaged 
resistance  may  be  input  if  na  is  set 
to  1.0.  Omit  Card(s)  14  if  IRES  a  1 

14 

1-8 

9-16 

• 

« 

JT 

Array  of  propeller  advanoe  coefficients 
(J)  in  asoending  order.  Maximum  of  60. 

10  points  per  oard. 

Omit  Card(s)  15  if  IOW  =  1 

15 

1-8 

9t16 

• 

KTO 

Array  of  propeller  thrusts  coefficients 
(K-j)  in  open-water  corresponding  to 
input  J'a.  Omit  Card(s)  16  if  IOW  =  1 

16 

1-8 

9-16 

* 

• 

KQO 

Array  of  propeller  torque  coefficient 
(Kq)  in  open-water  corresponding  to 
input  J's.  Omit  Card(s)  17  if  IOW  a  1 

17 

1-8 

9-16 

* 

NGP, 

Number  of  gear  ratios  considered 
--maximum  of  A 

18 

1-4 

NPROPS 

Number  of  sets  of  propellers  considers! 

18 

5-8 

GR 

Array  of  gear  ratios  (ra  ) 

If  m  -  0.0  is  input,  program 
will8compute  optimum  value  of  m 

19 

1-8 

9-16 

15 
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PROGRAM  PHPRLM 


DIN 

Propeller  diameter  (D)  in  inches 

Card 

20 

Columns 

1-8 

PD 

Propeller  pitch/diameter  ratio  (P/D) 

20 

9-16 

EAR 

Propeller  expanded  area  ratio  (EAR) 

20 

17-24 

Z 

Number  of  blades  per  propeller 

20 

25-32 

PEDES 

Brake  power  of  each  prime  mover 
used  for  sizing  propeller 

20 

33-40 

XI 

Index  of  ship  speed,  from  Card  7, 
for  sizing  propeller 

20 

41-48 

XJ 

Index  of  wave  height,  from  Card  8, 
for  sizing  propeller  and  for 
calculating  optimum  mg 

20 

49-56 

If  DIN  >  0,  do  not  input  PEDES  and  XI 

If  DIN  ■  0,  program  will  calculate  a 
diameter  based  on  above  speed,  wave 
height,  and  power 

Values  of  P/D,  EAR,  and  Z  must  be  input, 


Note;  Repeat  Card  20  for  each  set  of  propellers  considered. 
Then  repeat  Cards  2  through  20  for  each  case, 
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OUTPUT : 

Page  1  —  Echo  of  Input  Data 

Page  2  —  Characteristics  of  Prime  Movers  (Not  Printed  if  IPRT  > 0) 


MAX.BHP 


MAX.RPM 


GEAR  RATIO  mfl 


BHP/SHP 


SHP/DHP 


BHP  PER 
ENGINE 


Maximum  brake  horsepower  of  each  prime 
mover,  input  from  Card  4 

Maximum  rpm  of  prime  mover,  from  Card  4 

rpm  of  prime  mover/propeller  rpm 
from  Card  4 

gear  loss  ratio,  input  from  Card  4 

shaft  loss  ratio,  input  from  Card  4 

number  of  prime  movers  =  number  of 
propellers  input  from  Card  4 


points  from  engine  oharaoteristio  ourve, 
input  from  Cards  9  and  10 


ENGINE  RPM  Ne' 
TOTAL  DHP  PD' 


PROP.  RPM  N' 
Q-FT.LB  Q' 


total  delivered  horsepower  at  propellers 
Pe  npr  /K1  /K2 

propeller  rpm  ■  Ne/  mg 

Maximum  developed  torque  in  ft-lb 
33000  PD'/  (2  it  N') 


Page  3  —  Propulsion  Coefficients  and  Resistance 


LOA-FT 


LP-FT 


BPX-FT 


HT-FT 


DISPL-LB 
BETA -DEO 


a  overall  ship  length  in  ft,  input  from  Card  3 
a  projeoted  chine  length  in  ft,  from  Card  3 
a  maximum  beam  over  ohines  in  ft,  from  Card  3 
a  draft  at  transom  in  ft,  from  Card  3 
a  displacement  at  rest  in  lb,  from  Card  3 
a  deadrise  at  midships  in  deg,  from  Card  3 


IP  ><  v  -  -  V. 
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LCG-FT 

AG 

s 

CG  from  transom  in  ft,  from  Card  3 

VCG-FT 

KG 

= 

CG  from  baseline  in  ft,  from  Card  3 

XACC-FT 

Xaoo 

s 

distance  from  transom  in  ft,  at  whioh 
accelerations  are  oomputed,  from  Card  3 

C-LOAD 

ca 

— 

beam  loading  coefficient 

= 

A/<  P  g  Bpx3  )  =  V/Bpx3 

LP/V13 

Lp/V1/3 

3 

slenderness  ratio 

LP/BPX 

L/B 

B 

length-beam  ratio,  based  on  projected  chine 
length  and  maximum  chine  beam 

D-IN 

fin 

r 

propeller  diameter  in  inches,  from  Card  4 

P/D 

P/D 

= 

propeller  pitch  ratio,  from  Card  4 

BAR 

EAR 

= 

propeller  expanded  area  ratio,  from  Card  4 

NPR 

npr 

3 

number  of  propellers,  from  Card  4 

V-KT 

vK 

3 

ship  speed  in  knots,  input  from  Card  7 

V-LOA 

VkVl 

3 

speed-length  ratio  based  on  overall  length 

FNV 

Fny 

3 

speed-disDlaoement  ooeffioient 

V  /  (gV  1'3)  1/2 

SIGMA 

0 

3 

oavitation  number 

1  -T 

1-t 

= 

thrust  deduction  factor 

1-W 

1-w 

3 

thrust  wake  factor  =  torque  wake  faotor 

EA 

= 

appendage  drag  faotor 

1-t,  1-w, 

,  n 

a  generated  from  Subroutine  PRCOEF  if  IPC  nl 

1-t,  1-w, 

,  n 

a  input  from  Cards  11,  12,  13  if  IPC  a  2 

RB/W 

Hb/W 

3 

resistanoe-weight  ratio  for  bare  hull  in 
still  water 

Rb 

generated  from  Subroutine  PHRES  if  IRES  s  1 

*b 

input  from  Card  14  if  IRES  3  2 

W 

3 

weight  of  oraft  =  displacement  at  rest  A 

RA/W 

Ra/W 

= 

resistance-weight  ratio  for  appendaged  hull 
in  still  water 

= 

(Rb/w)  /  na 

18 
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program  phprlm 


SHPB  EHPh  =  effective  horsepower  of  bare  hull 

;  Rb  V  /  550. 


ge  4  — 
IPROP 


D-IN 

D-FT 

P/D 

EAR 

B LA  DEC 


EHPa  =  effective  horsepower  of  appendaged  hull 
■  Rq  V  /  550 

propeller  Open-Water  Characteristics  (not  printed  if  IPRT>0) 


Indicator  for  propeller  type,  from  Card  6 

- 

1  for  Qawn-Burrill  type 

s 

2  for  Newton-Rader  type 

S 

3  for  Wageningen  B-sorew  type 

Din 

= 

Propeller  diameter  in  inohes,  from  Card  20 

D 

s 

Propeller  diameter  in  ft  -  Din/12 

P/D 

s 

Propeller  pitch  ratio,  from  Card  20 

E.A.R 

a 

Propeller  expanded  area  ratio,  from  Card  20 

Z 

- 

Number  of  blades  per  propeller,  from  Card  20 

DEPTH-FT 

h0 

S 

depth  of  oenter  of  propeller  below 
waterline  in  ft  =  TA  +  0.75  D 

SIGMA/VSQ 

o/V2 

3 

constant  for  cavitation  number 

A 

s 

<PA  +  PH  -  Py) / (p/2) 

PA 

s 

atmospheric  pressure  =  2116  lb/ft^ 

vapor  pressure  =  36  lb/ft^ 

PV 

s 

PH 

= 

static  pressure  *  P  8  ti6 

AP-SQFT 

Ap 

s 

Projected  area  of  propeller  in  sq.  ft 
(  TT  D2/14 )  (EAR)  (1.067  -  0.229  P/D) 

vA 

S 

speed  of  advance  =  (1-w)  V 

V 

s 

ship  speed  in  ft/sec 

J 

J 

a 

propeller  advanoe  coefficient  *  VA/(nD) 

n 

a 

propeller  revolutions  per  second  (rps) 

KT 

Kt 

a 

propeller  thrust  coefficient  =  T/ (p  n2  D^) 

T  =  thrust  per  cropeller  in  1 
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PROGRAM  PHPRLM 


KQ 

Kq 

2  5 

propeller  torque  coefficient  =  Q/(p  n  D  ) 

Q 

torque  per  propeller  in  ft-lb 

Tables  of  open-water  ftp  and  Kq  values  generated  from 
Subroutine  OWKTQ  if  IOW  s  1  or  input  from  Cards  15-17  if 
IOW  =  2. 

EP 

no  : 

propeller  open-water  efficiency 
<J/2ir)  (Kj/Kq) 

KT/J2 

Kt/J2 

2  ? 

propeller  thrust  loading  ■  T/(p  D  ) 

KQ/J2 

kq/j2  - 

3  2 

propeller  torque  loading  ■  Q/(p  D  ) 

KQ/J3 

kq/j3  - 

2  3 

propeller  power  loading  ■  Qn/ (p  D  ) 

TC 

Tc 

■ 

2 

thrust  load  coefficient  ■  T/(Ap  Vq  p/2) 

Kf/r%  (Ap/D2)  (J2  +  4.84)] 

QC 

Qc 

■ 

torque  load  coefficient  -  Q./(Ap  VQ^R  p/2) 
KQ/lh  (Ap/D2)  (J2  +  4.84)] 

S.  7/S 

a0.7R/O  “ 

cavitation  number  baaed  on  velocity  at  0.7R 
cavitation  number  based  on  advance  velocity 

- 

J2/(J2  +  4.84) 

EPMAX 

no 

max 

maximum  propeller  efficiency 

Page  5  - 

Propeller  Cavitation  Characteristics  (Not  Printed  if  IPRT>0) 

IPROP, 

D-IN,  etc. 

Same  as  Page  4 

Tl,  T2, 

Ql,  etc. 

Constants  from  Subroutine  CAVKTQ 

SIGMA 

a  = 

cavitation  number  „ 

<PA  +  pH  -  Py)  /  (1/2  p  V‘  ) 

J 

J 

advance  coefficient  ■  VA/(nD) 

KT 

*T 

2  4 

thrust  coefficient  ■  T/(^n  D  ) 

KQ 

kq 

2  5 

torque  coefficient  »  Q/(^  n  D  ) 

Tables  of  IL 
from  Subroul 

and  K-  at  various  v'b  generated 
line  CAVKTQ 

KrtWfjrs,-  j/furwjirjMWlW 
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Page  6  -  Propeller  Sizing  (Not  Printed  if  Diameter  la  Input) 


VA  (FPS) 

VA 

m 

speed  of  advance  in  ft /sec 

V  (1-w)  where  V  1b  design  ship  speed 

P/D 

P/D 

propeller  pitch  ratio 

EAR 

EAR  " 

propeller  expanded  area  ratio 

NPR 

n  ■ 

pr 

number  of  propellers  ■  number  of  prime  movers 

BLADES 

Z 

number  of  blades  per  propeller 

DIN 

D . 
in 

propeller  diameter  in  inches 

MIN.  DIAM 

D  .  ■ 

min 

minimum  diameter  based  on  101  back  cavitation 
criteria  from  Gawn-Burrlll  propeller  series 

MAX.  DIAM 

D 

max 

maximum  diameter  based  on  maximum  propeller 
open-water  efficiency 

OPT.  DIAM 

D  _ 
opt 

optimum  diameter  selected 

DFT 

D 

propeller  diameter  in  ft 

SIGMA 

a  ■ 

cavitation  number  based  on  V. 

A 

KT/JSQ 

lyj2 

thrust  loading  per  propeller 

JT 

JT  ■ 

2 

advance  coefficient  at  K^,/ J 

KT 

kt 

thrust  coefficient  at 

KQ 

kq 

torque  coefficient  at 

EP 

no 

propeller  efficiency  at  JT 

*  after  Hq 

indicates  operation  at  J  above  peak  efficiency 

PC 

nD 

propulsive  coefficient 

T-LB 

T 

total  thrust  requirement  in  lb  at  design  speed 
and  wave  height 

C  after  T 

indicates  more  than  10%  back  cavitation 

*  after  T 

indicates  maximum  thrust  due  to  cavitation 

Q-FT.LB 

Q 

total  torque  in  ft-lb 

RPM 

N 

propeller  rpm 

DHP 

PD 

total  power  developed  at  propellers 

.  -i.  -b,»  >  **--  '  .  i>t  »v  -M"  » 
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Page  7  -  Powering  Requirements  and  Accelerations 


LOA-FT, 

LP-FT,  etc 

Same  as  Page  3 

V-KT 

vK 

* 

Ship  speed  in  knots,  input  from  Card  7 

H13-FT 

Hl/3 

: 

significant  wave  height  in  ft,  from  Card  8 

RW/W 

Rt/W 

= 

total  resistance-weight  ratio  in  seaway 
(Ra/W)  ♦  (Raw/W) 

Raw/W 

s 

added  resistance-weight  ratio  in  waves 

1.3  *  (Lp/V  ^3)  2’  Pfly  (Reference 

T 

s 

total  thrust  requirements  in  lb 

Rt/(1-t) 

KT/JSQ 

KT/J2 

= 

thrust  loading  per  propeller 

T/(nprpVA2  D2) 

JT 

Jip 

- 

propeller  advance  coefficient  corresponding 
to  KT/J2 

KT 

kt 

s 

propeller  thrust  coefficient  at  Jt 

KQ 

Kq 

= 

propeller  torque  coefficient  at  Jt 

EP 

^0 

- 

propeller  efficiency  at  JT 

<*  after 

indicates  operation  beyond  peak  efficiency, i.e. 
high  J 

Jt,  Kti  Kq,  Hq  generated  from 

Subroutine  PRINTP  as  funotion  of  Kt/J2 
and  o 

n 

= 

propeller  rps  =  VA  /  (JT  D) 

.  -LB 

T 

a 

total  thrust  in  lb  =  Kt  P  n2  D**  npr 
oheok;  T  =  Rt  /  (1-t) 

C  after  T  indicates  more  than  10  percent  back  cavitation 

*  after  T  indicates  that  thrust  required  exceeds 

maximum  thrust  limit  of  propeller  due  to 
oavitation  —  unless  the  propeller  can 
operate  in  the  fully  cavitating  range,  i.e., 
high  rpm,  low  Jt, 

Q-FT.LB  Q  =  total  torque  in  ft-lb  =  Kq  Pn2  npr 

#  after  Q  indicates  that  torque  required  exoeeds 

torque  limit  of  the  prime  movers  at  the 
required  rpm 


PROGRAM  PHPRLM 


RPM 

N 

s 

propeller  rpm  a  60  n 

*  after  N 

indicates  that  rpm  required  exoeeds  rpm 
limit  of  the  prime  movers 

DHP 

PD 

s 

total  horsepower  developed  at  propellers 

Q  n  /  550 

pE 

s 

total  effeotive  horsepower 

Rt  V  /  550 

PC 

= 

propulsive  coefficient  a  Pe/P£> 
check:  nD  =  n0  nH  nR 

nH 

= 

hull  effioienoy  a  (1-t)  /  (1-w) 

= 

relative  rotative  effioienoy  a  1.0  since 
torque  wake  assumed  equal  to  thrust  wake 

TRIM 

T 

trim  angle  in  deg 

fixed  trim  angle,  if  input  from  Card  3 
otherwise  trim  generated  by  Subroutine  SAVIT 

CG  ACC 

fiCG 

m 

average  1/10  highest  vertioal  accelerations 
at  the  oenter  of  gravity  in  g's 

7.0  (Hl/3/Bpx)  (1  +  t/2)0,25  (Lp/V173)"1,25  FnV 

BOW  ACC 

aBOW 

m 

average  1/10  highest  vertical  acceleration 
at  90£  of  Lqa  forward  of  transom  in  g's 

10.5  <H1/3/Bpx)  (1  +  t/2)0’5  W^V-75 

X  ACC 

ax 

m 

average  1/10  highest  vertical  accelerations 

at  location  X  in  g's 

acc 

■ 

*CG  +  <*»0H  -  *C0>  <*.«  -  LC0>'<0-9  l0A  -  LC0> 

Equations  for  aCG  and  aR0W  from  Reference  8. 


PROGRAM  PHPRLM 


;•  8  -  10  Percent  Back  Cavitation  (Gawn-Burrill  Propellers) 

V-KT  Vj,  -  ship  speed  in  knots,  from  Card  7 

SIGMA  -  cavitation  number 

T-LB  T  *  total  thrust  in  lb  at  which  10%  back 

cavitation  occurs 

rc  -  0.494  <rQt7R 

Q-FT.LB,  RPM,  DHP  -  corresponding  values  of  torque,  rpm^delivered 

power 

J,  KT,  KQ,  EP  ■  corresponding  propeller  characteristics 

;e  8  -  Maximum  Speed  with  Optimum  Gear  Ratio  (Input  m  ■  0) 

8 


V-KT 

VK 

■  ship  speed  in  knots  attainable  in  specified 

max 

sea  with  maximum  power 

H13-FT 

Hl/3 

■  significant  wave  height  in  ft  for  computing 
optimum  gear  ratio 

DHP 

PD 

*  maximum  power  delivered  at  propellers 

max 

RPM-P  N 


OPT. GEAR  m 
RATIO  8opt 


propeller  rpm  required  at  V„  and  H, ,, 


propeller  rpm  required  at  VR  and 

max 

optimum  gear  ratio 

engine  rpm  at  maximum  power 


PROGRAM  PHPRLM 


Page  9  —  Thrust  at  Maximum  RPM  of  Prime  Movers 


GEARRATIO  m 


SIGMA 


Q-FT.LB 


m  ■  gear  ratio ,  input  from  Card  19  or  optimum 

8  ra  computed  by  program 

8 

Nmax  =  maximum  rpm  of  propellers 

s  /  m„ 

®max  8 

nmax  =  Maximum  rps  of  propellers  a  Nmax/60 

Vj(  =  Ship  speed  in  knots,  input  from  Card  7 

o  =  oavitation  number 

Jmax  =  propeller  advance  coefficient  at  max.  rpm 
=  VA  /  (n^  D) 

Kt  *  propeller  thrust  coefficient  at  Jmax 

Kq  a  propeller  torque  ooeffioient  at  Jmax 

Hq  a  propeller  effiolenoy  at  Jmax 

Kji  Kq,  n0  generated  from  Subroutine 
PRINT?  as  function  of  Jmax  and  o 

Tmax  =  maximum  thrust  in  lb  available  at  Nmax 

■  «I  P"Ssx  O11  Ppr 

C  after  T  indicates  more  than  10%  back  cavitation 
max 

*after  T  indicates  limit  due  to  cavitation 
max 


Qmax  s  torque  in  ft-lb  at  Nmax 
s  kQ  P  tyRax  npr 

v after  U  indicates  that  torque  limit  of 
a  prime  movers  is  exceeded 

P_  »  power  delivered  at  propellers  • 


n _ /55C 


nax  max 


PROGRAM  PHPRLM 


Page  9  —  Thrust  at  Torque  Limits  of  Prime  Movers 


DHP 

PD' 

= 

developed  horsepower  1 

oaloulated  from  input 
points  for  engine 

RPM 

N' 

propeller  rpm 

characteristics —  See 
output  Page  2 

Q-FT.LB 

Q' 

n' 

- 

developed  torque 

propeller  rps  a  Nf/60 

V-KT 

V 

ship  speed  in  knots  at  whioh 

Kq  =  Q'/(  Pn*D5npr)  matches 

J’  3  V  7(n'D)  from  propeller  curves 
obtained  by  iteration  and  interpolation 

SIGMA 

o' 

= 

Cavitation  number  at  Vr' 

J 

J' 

s 

propeller  advanoe  ooeffioient  at  n' 

VA  ♦  /  <  n  •  D  > 

kt 

kt 

s 

propeller  thrust  ooeffioient  at  J' 

Kq 

Kq 

s 

propeller  torque  ooeffioient  at  J' 

EP 

no 

propeller  effioienoy 

%,  Kq,  rig  generated 
as  function  J'  and  d* 

at  J' 

front  Subroutine  PRINTP 

T-LB 

T' 

= 

thrust  in  lb  available=  KTp(n')2  npr 

Q-FT.LB 

Q' 

= 

torque  in  ft-lb  oheck  =  K^pCn')2  D5  npr 

DHP 

pd’ 

s 

developed  horsepower 

oheok  =  2 tt  Q*  n'/550 

Page  9  -- 

Speed-Power 

Limits 

H13-FT 

H 1/3 

= 

significant  wave  height  in  ft,  input  from 
Card  8 

V-KT 

VKmax 

" 

ship  speed  in  knots  at  whioh  thrust  required 
at  H^/3  (output  page  6 )  reaches  the  thrust 
limit  of  the  prime  movers  (output  page  7  ) 
due  to  either  rpm  or  torque  restriction, 

whichever  is  lower 


PROGRAM  PHPRLM 


T-LB 

T  * 

total  thrust  in  lb  at  V,. 

*Tnax 

C  after  T 

indicate*  more  than  10X  back  cavitation 

*  after  T 

indicates  thrust  limit  due  to  cavitation 

Q-FT.LB 

Q  » 

total  torque  in  ft-lb  at  V„ 

iSnax 

*  after  Q  indicates  torque  limit  of  primo  movers 

RPM 

N  s 

propeller  rpm  at  VR 

^max 

•  after  N 

indi oates  rpm  limit  of  prime  movers 

DHP 

PD 

total  developed  horsepower  at  Vv 

Nnax 

Page  10  — 

Habitability  Limits 

V-KT 

If 

ship  speed  in  knots,  input  from  Card  7 

1.0G 

l.Og  s 

average  1/10  highest  vertioal  accelerations 
representing  the  enduranoe  limit  of  the  orew 
for  a  maximum  of  4  to  8  hours 

1.5G 

1.5g  * 

average  1/10  highest  vertioal  aooelorations 
representing  the  enduranoe  limit  of  the  orew 
for  1  to  2  hours 

H 1 3 -FT 

«l/3 

signifioant  wave  height  in  ft.  corresponding 
to  habitability  limit  of  1.0  g  (or  1.5g)  — 
interpolated  from  aooelerations  at  location 

Xaoo  on  output  Page  6 

Note:  acceleration  predioatora  are 
considered  not  valid  for  H^3  >0.75  Bpx 
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NAME: 

SUBROUTINE  PRCHAR 

PURPOSE: 

Generate  propeller  open-water  and  cavitation 
characteristics.  Select  suitable  propeller 
diameter,  if  not  input 

CALLING  SEQUENCE: 

CALL  PRCHAR 

SUBPROGRAMS 

CALLED: 

OWKTQ,  CAVKTQ,  PRINTP,  TINT* 

INPUT: 

DIN 

Din 

propeller  diameter  in  Inches,  from  Card  20 

If  not  input,  a  suitable  diameter  will  be 
selected  by  this  routine. 

PD 

P/D 

■ 

propeller  pitch  ratio,  from  Card  20 

EAR 

EAR 

SB 

expanded  area  ratio,  from  Card  20 

Z 

Z 

m 

number  of  blades,  from  Card  20 

VADES 

VA 

a 

design  speed  of  advance  in  ft/ sec  derived 
from  design  ship  speed,  indexed  from  Card  20 

tdes 

T 

m 

total  thrust  in  lb  required  at  design  speed 
and  wave  height,  indexed  from  Card  20 

PDES 

PD 

at 

total  design  power  at  propellers  derived  from 
design  brake  power  input  on  Card  20 

PRN 

- 

number  of  propellers,  from  Card  6 

PROPELLER  OPEN-WATER  CHARACTERISTICS:  See  Subroutine  OWKTQ 
PROPELLER  CAVITATION  CHARACTERISTICS:  See  Subroutine  CAVKTQ 
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SELECTION  OF  PROPELLER  DIAMETER!  (When  D  is  not  input) 


DMIN 


DMAX 


DIAM 

DHP 

RPM 

DIN 


^min 


D 


max 


minimum  diameter  in  inches  based  on  10%  back 
cavitation  criteria  for  adequate  blade  area 
of  Gawn-Burrill  type  propellers 

“  V  1  >  C"„r  f  ^  (KI/j2,] 

,  _  0«88 
where  K^/J2  is  derived  from  tc  -  0.A9A 

representing  the  Gawn-Burrill  10%  cavitation 
line 


maximum  diameter  in  inches  based  on  maximum 
propeller  efficiency 

2 

same  equation  as  above,  with  K^/J  at  point 
of  maximum  rig 


array  of  diameters  in  increments  of  5  inches 


ft°"  t0  »■« 


array  of  power  requirements  for  above  diameters 


array  of  propeller  rpm  requirements  for 
above  diameters 


opt 


optimum  diameter  selected,  rounded  up  to 
next  full  inch 


D  .  if  power  required  <  design  power 
min 


is  interpolated  from  D  array  if  design  power  is 
between  requirements  for  D^^  and  Dmju( 


D  ^  •»  D  or  D  .  ,  whichever  is  larger,  if  both 

opt  mftx  min 

power  requirements  exceed  design  power 
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NAME: 


PURPOSE: 


SUBROUTINE  OWKTQ 

Calculate  propeller  open-water  characteristics  as 
function  of  pitch  ratio,  expanded  area  ratio,  and  number 
of  blades  from  coefficients  derived  from  Wageningen 
B-Scrc.w  Series. 


REFERENCE: 


Oosterveld  and  Van  Oossanan,  "Recent  Development  in 
Marine  Propeller  Hydrodynamics,"  Proceedings  of  the 
Netherlands  Ship  Model  Basin  40th  Anniversary  (1972), 
and  "Further  Computer  Analyzed  Data  of  the  Wageningen 
B-Screw  Series",  International  Shipbuilding  Progress, 
Vol.  22  (July  1975). 


CALLING  SEQUENCE:  .  CALL  OWKTQ 


INPUT: 


OUTPUT : 


■  propeller  pitch/diameter  ratio 

■  propeller  expanded  area  ratio 

■  number  of  propeller  blades 


'0 


flj.  -  number  of  J  values  generated  ~  max  of  60 

J  ■  array  of  propeller  advance  coefficients  in 

ascending  order  from  (J-0.)  to  (J  at  K»0.) 
in  increments  of  0.025  if  P/D£l.2 
in  increments  of  0.05  if  P/D^1.2 

K_  ■  array  of  open -water  thrust  coefficients 

-  f  (P/D,  EAR,  Z,  J) 

K  ■  array  of  open-water  torque  coefficients 

q  -  f  (P/D,  EAR,  Z,  J) 

Kt  and  K  developed  from  equations  in  above  references. 
For  Gawn^Burrill  type  propellers  (IPROP-1)  the  equations 
are  modified  to  produce  slightly  higher  anl  KQ 
than  the  Wageningen  B-Screw  Series,  ^ 


-,V' 


..  I  ■ij.-'W-  •'**>  i" 1 ■ 
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NAME:  SUBROUTINE  CAVKTQ 

TURPOSE:  Calculate  propeller  oharaoteriatioa  in  cavitation 

regime  as  function  of  pitoh  ratio,  expanded  area 
ratio  and  cavitation  number.  Generate  CALCOMP  plots 
of  Kt  and  Kq  versus  J. 

REFERENCE:  Blount  and  Fox,  "Design  Considerations  for 

Propellers  in  a  Cavitating  Environment,"  Marine 
Technology  (Apr  1973) 

CALLING  SEQUENCE:  CALL  CAVKTQ 


SUBPROGRAMS  CALLED : TQMAX ,  CALCOMP  Routines 
INPUT: 

J.PROP  Control  for  type  of  propellers 

=  1  for  Gawn-Burrill  type 

(flat  face,  segmental  sections) 

=  2  for  Newton-Rader  types 
a  3  for  Wageningen  B-Sorew  (airfoil  sections) 


PD 


P/D  =  propeller  pitoh/diameter  ratio 


EAR 

NJ 

JT 

KTO 

KQO 

NS 

SIOMA 
IP  LOT 


EAR  a  propeller  expanded  area  ratio 

nj  s  number  of  J  values  input  from  open-water 

curves  —  max.  of  60 

J  a  array  of  propeller  advance  coefficients 

K?  a  corresponding  array  of  propeller 

open-water  thrust  coefficients 

Kq  s  corresponding  array  of  propeller 

open-water  torque  ooeffioients 

na  a  number  of  oavitation  numbers  —  max.  of  8 

—  at  which  propeller  oharaoteristios  are 
to  be  oomputed  and  printed  from  this 
routine  (if  ns  s  0  only  the  constants 
are  computed) 

a  s  array  of  oavitation  numbers 

Control  of  CALCOMP  plots 

=  2  for  plots  of  Kj  and  Kq  vs  J  at  eaoh  a 
(no  plots  done  if  IPLOT  <  2) 


31 


SUBROUTINE  CAVKTQ 


GENERAL  NOTATION  FOR  PROPELLERS: 

Va  =  propeller  speed  of  advance 

n  =  rate  of  revolution 

D  =  propeller  diameter 

T  =  thrust 

Q  s  torque 

p  =  water  density 

P0  =  pressure  at  center  of  propeller  s  pa+PH"PV 

J  s  advance  coefficient  =  Va/  (n  D) 

Kf  =  thrust  ooeffioient  =  T  /  (  p  n2  D^) 

Kq  s  torque  ooeffioient  s  Q  /  ( p  n2  D5) 

K-p/J2  =  thrust  loading  =  T  /  (  P  D2Va2) 

Kq/J2  =  torque  loading  -  Q/(P  D3  Va2^ 

Kq/j3  s  power  loading  a  Q  n/  (  p D2  Va^) 

o  =  oavitation  number  based  on  advanoe  velooity 

=  P0  /  (1/2  pVA  2) 

V0  7R2  =  velocity  2  at  0.7  radius  of  propeller 

=  Va2  +  (0.7  n  nD)2  =  Va2(J2+4.8'I)/.I2 

0o,7R  =  cavitation  number  based  on  Vq,7R 

=  P0/(l/2  PV0-7R2)  =  0 j2/(j2+4.84) 

Ap  =  projected  area  of  propeller 

=  (  tt  D2/4)  EAR  (1.067-0.229  P/D) 

s  thrust  load  ooeffioient 
=  T  /  C 1/2  p  Ap  Vo,7R  2) 

=  KT  /  1/2  (Ap/D2)  (J2+H.8iO| 

Q0  =  torque  load  coefficient 

=  Q  /(1/2  p  Ap  Vq, 72R  2) 

=  Kq  /  [l/2  (Ap/D2)  (J2+1*.8H)1 


SUBROUTINE  CAVKTQ 


MAXIMUM  THRUST  AND  TORQUE  LOADS: 


OUTPUT: 

T1 


Blount  and  Fox  (see  referenoe)  give  equations  for 
maximum  thrust  and  torque  load  ooeffioients  in  a 
oavitating  environment  based  on  regression  of 
experimental  data  for  the  three  propeller  series  used 


herein. 

T°m 

= 

maximum  thrust  load  coefficient 
a  0Q,7Rb  (transition  region) 

= 

Tox( fully  oavitating  region) 

Q°m 

s 

maximum  torque  load  ooeffioient 
o  (transition  region) 

^ox  (fully  oavitating  region) 

IPROP 

a 

a 

1.2 

1 

a 

a 

0.703  +  0.25  P/D 

2 

a 

a 

1.27 

3 

I- 

iti 

I 


T2 

b 

5 

1.0 

1 

b 

5 

0.65  +  0.1  P/D 

2 

b 

: 

1.0 

3 

Q1 

0 

s 

0.200  P/D 

1 

c 

s 

0.2U0  P/D  -  0.12 

2 

0 

s 

0.247  P/D  -  0.0167 

3 

Q2 

d 

3 

0.70  +  0.31  EAR0 

1 

d 

3 

0.50  +  0.165  P/D 

? 

d 

3 

1.04 

3 

TCX 

T°x 

S 

0.0725  P/D  -  0.0340  EAR 

1 

T  * 

s 

0.0833  P/D  -  0.0142  EAR 

2 

T*x 

3 

0.0 

3 

QCX 

3 

[0.0185  (P/D)?  -  0.0166  P/D  + 

/EAR^J 

0.00594] 

<3ox 

3 

0.0335  P/D  -  0.024  EAR1/2 

2 

y«x 

3 

0.0 

3 

RMAX 

k 

a 

0.8 

Sinoe 

full- 

•scale  trial  data  (see  Figures 

5  and  6  of 

referenoe)  indicates  actual  thrust  and  torque  in  the 
transition  region  less  than  the  maximums  derived  from 
the  propeller  series  data,  the  faotor  k  is  applied  to 
TCm  and  QCm  in  the  transition  region.  The  faotor 
k  is  not  applied  to  T0x  and  Q0x. 


33 

l 


f 


SUBROUTINE  CAVKTQ 

APD2 

Ap/B2/2 

= 

Constant  for  calculation  of  T0  and  Q0 

J 

J 

= 

advance  coefficient  from  input  array 

OPEN  WATER  ] 

- 

input  values  of  open-water 

KT  KQ  j 

thrust  and  torque  ooeffioients 

SIGMA 

a 

= 

cavitation  number  from  input  array 

KT 

Kp 

£ 

thrust  coefficient  as  f  (J,  o ) 

3 

Kt  or  Kt  ,  whichever  is  smaller 
*o 

KTm 

= 

T0n  (1/2  Ap/D2)  (j2  +4.84) 

T°m 

S 

(k  a  o0#7R»>)  or  (t0x), 

whichever  is  greater 

LC 

. 

1  character  identifier  for  propeller 

cavitation 

C  indicates  more  than  10X  back  cavitation 
for  Gawn  propst  tc>  0.494  Oq 

*  indicates  thrust  limit  due  to  cavitation 


KQ 


m 


=  torque  ooeffioient  as  f  (J,  o  ) 

=  Ko  or  Kn  i  whichever  is  smaller 
'■'0  wm 


Qc 


(1/2  AD/D2)  (j2  +  4.84) 

m  r 


(k  o  o0,7Rd)  or  (Q0x) * 
whichever  is  greater 

Kxm  and  KQffl  generated  by  Function  TQMAX 


CALCOMP  PLOT? :  If  IPL0T=2,  open-water  Ky  and  Kq  as  well  as  ftp  and 
Kq  representing  the  transition  region  and  fully 
oavitating  region  at  eaohaare  plotted  as  a  function 
of  J. 
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NAME: 
PURPOSE : 

CALLING 

INPUT: 

SIGMA 

JT 

ITQ 


OUTPUT: 

TQMAX 


FUNCTION  TQMAX 

Calculate  maximum  thrust  or  torque  coefficient  in  a 
oavitating  environment  as  funotion  of  oavitation 
number  and  advance  coefficient 

SEQUENCE:  X  =  TQMAX  (SIGMA,  JT,  ITQ) 


o  s  oavitation  number 

J  =  advance  coefficient 

i  =  1  if  maximum  thrust  coefficient  required 

1  =  2  if  maximum  torque  oooffioient  required 


Variables : 


2 

a,  b,  c,  d,  tCj{,  QCx,  k,  1/2  Ap/D 
generated  by  Subroutine  CAVKTQ 


or  K„ 


: 

s 


depending  on  value  of  i 


maximum  thrust  load  ooeffioient 
k  a  O0>7.Rb,  or  Tc.x  if  greater 

cc  (1/2  Ap/D2)  (J2+1J.810 
m 


Q0jn  =  maximum  torque  load  ooeffioient 
S  k  C  °0.7R  *  or  if  grater 

KQ  =  QCm  (1/2  Ap/D 2)  (J‘%4.84) 

v'tn 
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NAME: 


SUBROUTINE  PRINTP 


PURPOSE:  Interpolate  for  propeller  performance  at  specified 

value  of  (1)  advance  coefficient  J,  (2)  thrust 
loading  Kp/J2,  (3)  torqut  loading,  Kq/J2,  or 
(4)  power  loading  Kq/j3. 

CALLING  SEQUENCE:  CALL  PRINTP  (IP,  PCOEF,  SIGMA) 


SUBPROGRAMS: 

TQMAX, 

YINTE 

INPUT: 

IP 

Option 

*  1* 

2,  3,  or  4 

PCOEF 

= 

input  propeller  coefficient, 
dependent  on  value  of  IP 

jt 

s 

advance  coefficient,  input  if  IP=1 

KT/j2 

s 

thrust  loading,  input  if  IPs 2 

Kq/J2 

s 

torque  loading,  input  if  IPs 3 

Kq/j3 

s 

power  loading,  input  if  IPs4 

SIGMA 

a 

s 

oavitation  number 

NJ 

nJ 

a 

number  of  J  values  defining  propeller 
characteristics 

JT 

J 

a 

array  of  advance  ooeffioient,  in 
ascending  order 

KT 

kt0 

= 

array  of  open-water  thurst 
coefficients 

KQ 

PERFORMANCE  AT 

K«o 

SPECIFIC  J: 

array  of  open-water  torque 
coefficients 

JTP 

Jp 

= 

input  advanoe  ooeffioient 

KTr 

Kt 

a 

thrust  ooeffioient  at  Jp 
open-water  thrust  ooeffioient 
interpolated  from  input  array  of 

K>p  versus  J,  or  maximum  thrust 
ooeffioient  in  oavitating  regime 

Kp  calculated  by  Funotion 

TQMAX,  whichever  is  smaller, 

KQP 

Kq 

= 

torque  ooeffioient  at  Jp 
open-water  value  interpolated  from 

Kq  vs  J,  or  maximum  oavitation 
value  Kq  calculated  from  TQMAX, 
whioheveP  is  smaller 
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SUBROUTINE  PRINTP 


PERFORMANCE  AT  SPECIFIC  LOADING: 


plog 

ln(KT/J2) 

if 

IP=  2  I 

natural  log  of 

ln(KQ/j2) 

if 

IPs  3 

input  loading 

ln(KQ/j3) 

if 

IPs  4  1 

coefficient 

XLOG 

ln(KT  / J?) 

ln(KQo/Jy 

ln(Kq°/j3) 

if 

IP=2 

array  of  natural  logs 

if 

IP=3 

of  open-water  loading 

if 

IPs  4  j 

ooeffioient  at  J  value 
from  input  array 

JTP  Jt  :  open-water  advance  coefficient  interpolated 

0  from  array  of  open-water  loading 

coefficients  versus  J  at  the  speoifio 
loading  required  (logs  are  used  because  of 
the  rapid  ohange  of  loading  ooefrioient  at 
low  J*s) 

If  Jt0  is  in  non-cavitating  region  (Kt0<  KTm) 

KTP  Kt  thrust  and  torque  coefficients  at  Jtq 

KQP  Kq  interplated  from  arrays  of  Ktq  and  Kqq  vb  J 


If  Jtq  is  in  oavitating  region  (Kt0  > 

XLOG  ln^KT_  /J^)  IP®2  I  array  of  natural  logs 

ln(Ky“J  /J2)  if  IP=3  J  of  loading  ooeffioients 

ln(KQ‘"  /j3)  if  ip-n  j  based  on  Kt  or 

Kq  as  function  J 

JTP  Jt  =  advanoe  coefficient  interpolated  from  array 

“  of  oavitation  loading  coefficients  vs  J  at 

the  speoifio  loading  required 

KTP  KT  \  maximum  oavitation  thrust  and  torque 

KQP  KQ  /  coefficients  at  Jt  calculated  from  TQMAX 

m 

OUTPUT: 

JTP  Jf  a  final  advanoe  coefficient  \ 

KTP  Kt  a  final  thrust  coefficient  I  at  propeller 

(  performance 

KQP  Kq  a  final  torque  coefficient  /  point 

I  specified  by 

EP  n0  =  propeller  efficiency  I  PCOEF  and 

=  JT  Kt/ (2 tt  Kq)  I  SIGMA 
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•  4* 


XSIG7  4  94  o  0,88 
0.7R 


SUBROUTINE  PRINTP 


thrust  load  coefficient 

Kj  /  Da(Vc2)  <j2+4‘84>] 

cavitation  number  based  on  velocity 
at  0.7  radius  of  propeller 

o  J2/  (J2+4,84)  4.84-(0.7t t)‘ 


■  term  representing  10Z  back  cavitation  line  for 
Gawn-Burrill  propeller  series 

■  1  character  identifier  for  propeller  cavitation 
*  indicates  thrust  limit  due  to  cavitation: 


C  indicates  more  than  10Z  back  cavitation  for 
Gawn-Burrill  propellers,  but  less  than 
thrust  limit  cavitation 
0.8B 

Tc  >  0.494  aojR 


-i'lWk  ' -tlWt  •  iil  ■  ’’fi  Jt  • 


NAME: 


SUBROUTINE  TPLOT 


PURPOSE: 


Generate  CALCOMP  plots  of  thrust  versus  speed  for 

(1)  thrust  required  in  various  sea  states 

(2)  thrust  limit  at  maximum  rpm  of  prime  movers 

(3)  thrust  at  torque  limit  of  prime  movers 
with  one  or  more  gear  ratios 

(4)  thrust  at  10%  back  cavitation,  if  IPROP  ■  1 


CALLING  SEQUENCE:  CALL  TPLOT 


SUBPROGRAMS 


CALCOMP  Routines,  FACTOR,  PLOT,  AXIS,  SYMBOL , 
NUMBER,  FLINE,  DASHL 


INPUT: 

TITLE 


■  identification  to  be  printed  at  top  of 
graph  80  characters  maximum 


■  array  of  ship  speeds  in  knots 

■  array  of  significant  wave  heights  in  ft. 

*>  number  of  speeds  —  maximum  of  20 

■  number  of  wave  heights  ~  maximum  of  5 

•  number  of  points  defining  torque  limits  of 
prime  movers  --  maximum  of  10 

■  number  of  gear  ratios  —  maximum  of  4 

■  matrix  of  values  for  thrust  required  as 
function  of  and 

-  array  of  thrust  values  representing  10% 
back  cavitaiton,  if  IPROP  -  1 


■  array  of  thrust  values  at  maximum  rpm  of 
prime  movers 

■  array  of  speeds  in  knots  at  which  torque 
limits  are  defined 

■  array  of  thrust  values  corresponding  to 
torque  limits  of  prime  movers 

*  array  of  gear  ratios 


rrsnr 


OUTPUT: 

X-axes 

Y-axes 

Carves 


SUBROUTINE  TPLOT 


Ship  speed  V  in  knots 
Ship  speed  V  In  m/sec 

Thrust  in  lb 
Thrust  in  N 

(1)  Thrust  required  at  each  wave  height 

(2)  Thrust  at  engine  rpm  limit,  for  each  gear  ratio 

(3)  Thrust  at  engine  torque  limit,  for  each  gear  ratio 

(4)  Thrust  representing  10%  back  cavitation  (dash  line) 


Intersection  of  Curve  (1)  and  the  lower  of  curves  (2) 
and  (3)  represents  the  maximum  speed  obtainable  at  each 
wave  height  and  gear  ratio. 
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NAME 


SUBROUTINE  HPLOT 


PURPOSE! 


SUBPROGRAMS : 


INPUT! 


•  •  •  •  V  .aatM 


TITLE 


Generate  CALCOMP  plots  of  significant  wave  height  versus 
spaed  corresponding  to 

(1)  endurance  limit  of  crew  for  4  to  8  hours  operation 

(2)  endurance  limit  of  crew  for  1  to  2  hours  operation 

(3)  power  limit  of  propulsion  system  for  1  or  more 

gear  ratios 


CALLING  SEQUENCE:  CALL  HPLOT 


CALCOMP  Routines  FACTOR,  PLOT,  AXIS,  SYMBOL,  NUMBER,  FLINR 


■  identification  to  he  printed  at  top  of  graph 
maximum  of  80  characters 

■  array  of  ship  speeds  in  knots 

■  array  of  significant  wave  heights  in  ft 

■  number  of  speeds  —  maximum  of  20 

*  number  of  wave  heights  —  maximum  of  5 

■  number  of  gear  ratios  —  maximum  of  4 

•  array  of  significant  wave  heights  as  function 
of  Vj.  corresponding  to  average  1/10  highest 
vertical  accelerations  of  l.Og 

■  array  of  significant  wave  heights  as  function 
of  V„  corresponding  to  vertical 
acceleration  of  1.5g 

■  array  of  maximum  ship  speeds  in  knots  as 
function  of  H,  derived  from  intersection 
of  required  thrust  at  each  wave  height  with 
thrust  curve  at  engine  rpm  or  torque  limit, 
whichever  is  lower 

■  distance  forward  of  transom  in  ft  at  which 
accelerations  have  been  computed 

■  array  of  gear  ratios 


imi, 


r 


OUTPUT i 

X-axes 

Y-axes 

Curve a 


SUBROUTINE  HPLOT 

Ship  speed  V  in  knots 
Ship  spaed  V  in  m/sec 

Significant  wave  height  H.  ,3  in  ft 
Significant  wave  height  in  m 

Wave  height  envelopes  representing  average 
1/10  heighest  vertical  accelerations  of 

(1)  1,0  g  —  habitability  limit  for  4  to  8  hours 

(2)  1.5  g  —  habitability  limit  for  1  to  2  hours 

(2)  1  Envelopes  of  maximum  speed  versus  wave  height 
(4)  >  due  to  power  limit  of  the  propulsion  system, 
etcj  for  1  to  4  gear  ratios 


NAME: 

PURPOSE: 

REFERENCE: 

CALLING  SEQUENCE! 

SUBPROGRAMS: 

INPUT: 

FNV 

OUTPUT : 

TDF 

ADF 

TWF 

PROCEDURE: 


SUBROUTINE  PRCOEF 

Estimate  propulsion  ooeffioients  for  planing  hull 
with  propellers  on  inolined  shafts 

Blount,  D.L.  and  D.L.  Fox,  "Small  Craft  Power 
Predictions,"  Western  Gulf  Seotion  of  the  Society  of 
Naval  Arohiteots  and  Marine  Engineers  (Feb  1975) 

CALL  PRCOEF  (FNV,  TDF,  ADF,  TWF) 

MINP,  YINTE 


Fn v  s  speed-displacement  coefficient 

«  V/(g  V 1 '3)1/2 


1-t  a  thrust  deduction  factor 

a  total  horizontal  reslstanoe  (Rt) 
total  shaft-line  thrust  (T) 

na  =  appendage  drag  factor 

a  resistanoe  of  bare  hull  (Rk) _ 

resistanoe  of  appendaged  hull  (Ra) 

1-w  a  thrust  wake  faotor  a  torque  wake  factor 

1-t,  1-w,  and  ha  interpolated  from  following  tabU  of 
values  at  input  value  of  Fn^  ,  The  tabulated  data 
represent  mean  values  from  a  bandwidth  of  data  aolleoted 
for  numerous  twin-sorew  planing  oraft  and  reported  in 
above  reference. 


^3 

<1 

1-t 

1-w 

Oa 

0.5 

0.92 

1.05 

0.951 

1.0 

0.92 

1.06 

0.948 

1.5 

0.92 

1.04 

0.942 

2.0 

0.92 

0.99 

0.934 

2.5 

0.92 

0.97 

0.925 

3.0 

C.92 

0.975 

0.913 

3.5 

0.92 

0.98 

0.900 

4.0 

0.92 

0.975 

0.885 

A3 


NAME: 


SUBROUTINE  PHRES 


PURPOSE:  Estimate  the  bare-hull,  smooth-water  resistance  of  a 

hard-ohine  planing  hull  from  synthesis  of  Series  62 
and  65  experimental  data 

REFERENCE:  Hubble,  E.N.,  "Resistance  of  Hard-Chine,  Stepless 

Planing  Craft  with  Systematic  Variation  of  Hull  Form, 
Longitudinal  Center  of  Gravity,  and  Loading,"  DTNSRDC 
Report  4307  (Apr  1974). 

Figures  9  and  10  of  DTNSRDC  Report  SPD-0840-01  (Deo 
1978,  Revised  Aug  1979) 

CALLING  SEQUENCE:  CALL  PHRES  (DLBS,  FNV,  SLR,  DCF,  SDF,  RH02,  VIS,  RLBS) 


SUBPROGRAMS: 

DISCOT, 

YINTX,  C1DSF 

INPUT: 

DLBS 

A 

* 

ship  displacement  in  lb 

FNV 

Fny 

8 

speed-displaoement  ooeffioient  (V/ (g/V1^3)1^2 

SLR 

Lp/  V 1^3  a  slenderness  ratio 

DCF 

aof 

»  correlation  allowance!  may  be  0.0 

SDF 

SDF 

a 

0.0  corresponds  to  mean  resistanoa-woight  R/W 
ourves  derived  from  Series  62  and  65  data 

- 

1.645  corresponds  to  minimum  R/W  ourves 

can  be  varied  to  approximate  R/W  for  a 
particular  hull  form  from  model  experiments 

RH02 

P/2 

- 

0.5  x  water  density  in  lb  x  seo^/rt^ 

VT.S 

V 

8 

water  viscosity  in  ftz/aeo 

OUTPUT: 

RLCS 

Rb 

8 

bare-hull,  smooth-water  resistance  in  lb 

8 

(mean  R/W  -  SDF  x  0 )  x  A 

0 

8 

standard  deviation  of  Series  62-65  data  from 
mean  R/W 

PROCEDURE: 

XFNV  array 

Tabulat 

ed  values  of  Fny  from  0.0  to  4.0 

ZSLR  array  Tabulated  values  of  Lp/V*^  from  4.0  to  10. 0 


SUBROUTINE  PHRES 

YRWM  matrix  Tabulated  values  of  mean  R/W  as  f(Fny  >  Lp/V^) 

for  100, 000- lb  planing  oraft  derived  rrom  Series  62  and 
65  experimental  data.  See  Table  1. 

YWSR  matrix  Tabulated  values  of  mean  wetted  area  ooeffioients  S/V2^3 
from  Series  62  and  65  hulls.  See  Table  2. 

SD  array  Tabulated  values  of  standard  deviation  a  as  f(Fnn) 

See  Table  1. 

RWM  R/W  for  100,000-lb  planing  oraft  interpolated  from  YRWM 

matrix  of  mean  R/W  values  at  input  Fn^  and  Lp/V 1^3 

WSR  S/ V  2^3  interpolated  from  YWSR  matrix  at  input  Fno 

andLp/vl'3 

Subroutine  D1SC0T  used  for  the  double  interpolation 

SDM  a  interpolated  from  SD  array  at  input  Fny 

Function  YINTX  used  for  single  interpolation 

RWM  (R/W)m  =  oorreoted  R/W  for  100,000~lb  planing  oraft 

«  (mean  R/W  interpolated)  -  (SDF  x  a  ) 

DLBM  Am  o  displacement  of  100,000-lb  planing  oraft 

XL  X  s  linear  ratio  of  aotual  ship  to  100,000-lb 

a  <A/Am)^3 

VFPSM  Vm  =  speed  of  100,000-lb  oraft  in  ft/seo 

»  (input  Fny)  x  19*32 

VFPSS  V «j  =  speed  of  aotual  ship  in  ft/seo  =  Vm  X1^2 

PLM  Lm  s  length  of  100,000-lb  oraft  in  ft 

a  11.6011*  (input  Lp/  V 3) 

PLS  La  =  length  of  aotual  ship  in  ft  a  Lm  X 

REM  Rnm  b  Reynolds  number  of  iOO, 000-lb  oraft 

*  vm  W  vm 

RES  Rna  3  Reynolds  number  of  aotual  ship  s  Va  La/vs 

CFM  Cp  b  Sohoenherr  friotional  resistance  ooeffioient 

for  100 ,000- lb  oraft 
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i) 


■Hull* 


SUBROUTINE  PHRES 


CFS 

CFS 

Sohoenherr  friotional  resistance  coefficient 
for*  aotual  ship 

Function  C1DSF  used  to  obtain  Sohoenherr 
friotional  resiatanoe  coefficients 

SM 

sra 

S 

wetted  area  of  100,000-lb  oraft  in  ft2 
134.5925  S/  V2/3 

SS 

S3 

= 

wetted  area  of  aotual  ship  in  ft2  sSn  X ^ 

RM 

Rm 

3 

resistance  of  100,000-lb  oraft  in  lb 
(R/W)n  An 

CTM 

CTtn 

3 

3 

total  resistance  ooeffioient  of  100, 000-lb 
oraft 

*m  /<V  Sra  pn/2) 

CR 

Cr 

3 

residual  resiatanoe  ooeffioient  s 

CV  cFm 

CTS 

cTs 

3 

total  resistance  ooeffioient  of  aotual  ship 
CjTg  Cp  +  ACp 

VIS 

Vs 

3 

kinematio  visoosity  for  aotual  ship 

VISM 

Vn 

a 

kinematic  visoosity  for  tabulated  data  s 
1.2817  x  10“5 

RH02 

/•s'* 

3 

1/2  water  density  for  aotual  ship 

RH02M 

CM 

E 

9 

1/2  water  density  for  tabulated  data 

1.9905/2 

RLDS 

Rb 

3 

resiatanoe  of  aotual  ship  in  lb 

*  cTa  V82  Ss  Pa/2 
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NAME:  SUBROUTINE  SAVIT 

PURPOSE:  Estimate  the  bare-hull,  smooth-water  resistance  and 

trim  for  a  hard-chine  planing  hull  using,  Savitsky's 
equations  for  prismatio  planing  surfaces 

CALLING  SEQUENCE:  CALL  SAVIT  (DISPL,  LCG,  VCG,  VFPS,  BEAM,  BETA,  TANB, 

COSB,  SINB,  HW,  WDCST,  RHO,  VIS,  AG, 
DELCF,  R,  TD,  NT,  CLM,  GDB) 

SUBPROGRAM:  CIDSF 

INPUT: 


DISPL 

A 

= 

ship  displacement  in  lb 

LCG 

AG 

s 

distance  of  oenter  of  gravity  CG  forward 
of  transom  in  ft 

VCG 

KG 

= 

distanoe  of  CG  above  baseline  in  ft 

VFPS 

V 

= 

speed'  in  ft/seo 

BEAM 

B 

- 

beam  in  ft 

a 

maximum  ohlne  beam  Bpg 

BETA 

8 

a 

deadrise  angle  in  degrees 

a 

deadrise  at  midships  Bm 

TANB 

tan  8 

COSB 

QOS  8 

SINB 

sin  8 

HW 

Hw 

= 

hoight  of  oenter  of  wind  drag  above 
baseline  in  ft 

WDCST 

CDW' 

a 

horizontal  wind  foroe  in  lb  /V2 

RHO 

P 

a 

water  density  in  lb  x  seo2/ft^ 

VIS 

V 

= 

kinematio  visoosity  of  water  in  ft2/seo 

AG 

g 

a 

acceleration  of  gravity  in  ft/seo2 

DELCF 

ACjp 

- 

correlation  allowance;  may  be  0 

OUTPUT : 

R 

Rb 

s 

bare  hull,  smooth-water  resistance  in  lb 
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SUBROUTINE  SAVIT 


TD 

T 

= 

trim  angle  in  degrees 

NT 

Number 

of  iterations  to  obtain  trim  angle 

CLM 

A 

= 

mean  wetted  length-beam  ratio  LD/B 

GDB 

2S 

= 

longitudinal  canter  of  pressure,  distanoe 
forward  of  transom,  in  ft 

PROCEDURE : 

TD 

T 

= 

trim  angle  of  planing  surfaoe  from 

horizontal  in  deg 

first  approximation  oft  a  deg 

CV 

cv 

: 

speed  coefficient  =  V/ (gB) 1 

CLM 

A 

s 

mean  wetted  length-beam  ratio 

Lm/B  =  (LK  +  LC)/2B 

CIO 

C^o 

a 

lift  coefficient  for  flat  surface 

XiA(0.0U  AI/2  +  0.0055  A5/2/Cv2) 

CLB 

-  lift  coefficient  for  deadrise  surfaoe 
=  A/[V2  Bzp/2j=  CLo-0.0065  CLq 

CL0  and  obtained  by  Newton-Raphson  iteration 
first  approximations:  Cl0  =  0.085;  A*  1.5 

XK 

f'K 

= 

wetted  keel  length  in  ft 

B[A  +  tang  /2ff  tanx  )j 

xc 

LC 

S 

wetted  chine  length  in  ft  =  2BA  -Lj( 

Lk-  Lq  s  (Btan  g  )/  (  tt  tan  x  ) 

GDB 

AP 

: 

longitudinal  oenter  of  pressure  forward  of 
transom  in  ft 

B  A [0 .  7 5  -  ]/(5. 21  Cv  /  >2  +  2.39)] 

CLD 

CLd 

~ 

dynamic  component  of  lift  coefficient 

0.012  A1^  x1.1 

VM 

\ 

a 

mean  velooit.y  ever  planning  surfaoe  in  ft/seo 

* 

V  [1_^cLd  '  °-0065  0  CLd0,6)/(A  cost)]  1/2 

PE 

R 

n 

= 

Reynolds  number  for  planing  surface 

Vn  B  A  /  v 

SUBROUTINE  SAVIT 


CF 

C  +  AC 

Sohoenherr  frictional  rosia canoe  ooeffoient 

*•  r 

as  f(Rn)  plus  correction  allowance 

DFX 

Dp 

= 

visoous  foroe  due  to  wetted  surface, 
parallel  to  the  planing  surface,  in  lb 

= 

(Cf+ACf)  (P/2  (Vm  2)  (B2  A/oos g  ) 

CK 

r‘K 

= 

1.5708  (1  -  0.1788  tan2goos8  -0.096H6 
tang  sin2 6) 

CK1 

cKi 

a1 

= 

Cj(  tant  /sing 

A1 

s 

[Hin2T(l-2CK)  +  CK2tau2T(l/8in2g-8in2T)]1/2 

cos  T  +  CF  tan  T  sin  T 

TAN 

tun  0 

= 

(&1+Ck1  )/(l-ai  C^) 

THETA 

0 

s 

angle  between  outer  spray  edge  and  keel  in 
radians 

* 

= 

arotan(tan<f>  oos  g  ) 

DLH 

AA 

= 

effective  inorease  in  length-beam  ratio  due 
to  spray 

c 

[tang  (  it  tant)j  -  l/(  2  tane  )  / (2  oos0  ) 

RE 

Rnc 

5 

Reynolds  number  of  spray 

S 

S 

VB  /(3oosg  sin  6  )  /V 

CF 

% 

s 

Sohoenherr  frictional  resistance  coefficient 
for  spray  drag 

DSX 

Ds 

= 

viscous  foroe  due  to  spray  drag,  parallel  to 
the  planing  surfaoe,  in  lb 

= 

CFs  ( P/2 )  V2)  (B2AA  /  oos 6  ) 

DVJX 

Dw 

= 

component  of  wind  drag  parallel  to  planing 
surfaoe  in  lb 

a 

Cj^,*  v2/oost 

DTX 

D*p 

a 

total  drag  foroe  parallel  to  planing  surfaoe 

in  lb 

= 

Dp  +  Dg  +  Dy 

PDBX 

PT 

a 

total  pressure  foroe  perpendicular  to 
surfaoe  in  lb 

A/oos'c  +  Df  tarn 


SUBROUTINE  SAVIT 

®p  «  moment  arm  from  oenter  of  pressure  to  CO  in  ft 

«  SB  -  JiP 

fp  a  moment  arm  from  center  of  vlaooua  force  to 

CO  in  ft 

■  ES  -  (B  tan  B /A) 

fw  a  moment  arm  from  center  of  wind  drag  to  CQ 

in  ft 

a  K5  -  HW 

EM  c  sura  of  moments  about  CQ  in  ft-lb 

a  Pf  Op  +  (Dp+Dg)  fp  +  Dyfy 

Iterate  with  small  ohariges  in  T  until  EMS0.001A 

Number  of  iterations  required  to  obtain  equilibrium 
trim;  maximum  of  15  iterations 

R  a  net  horizontal  resistance  Toroe  in  lb 

*  D<t>  cost  +  P«r  slnx 


NAME: 


FUNCTION  C1DSF 


PURPOSE : 

CALLING 

INPUT: 

XN1RE 
OUTPUT : 
C1DSF 


Caioulate  Sohoenherr  frictional  resiatanoe  coefficient 
SEQUENCE:  CF  =  C1DSF  (XN1RE) 


Rn  =  Reynolds  number  =  V  L  /  v 


Cp  s  Sohoenherr  friotional  resiatanoe  ooeffioiant 

Iteration  with  Newton-Raphaon  method 
Sohoenherr  formula:  0.242  / 8  I0810  ^F 


PROCEDURE: 


NAME: 


FUNCTION  SIMPUN 


PURPOSE:  Numerical  integration  of  area  under  curve  defined  by 

set  of  (x,y)  points  at  either  equal  or  unequal 
intervals 

CALLING  SEQUENCE:  AREA  =  SIMPUN  (X,  Y,  N) 

INPUT : 

X  array  Table  of  x  values— independent  variable  - 

x  values  must  be  in  asoending  order 

Y  array  Table  of  y  values— dependent  variable 

N  Number  of  (x,y)  values 

OUTPUT: 

SIMPUN  Area  under  curve*  / y  dx 


.»h«  iiiituMwini11  IfWTOHBli'ig 
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NAME: 


PURPOSE: 


CALLING  SEQUENCE: 


SUBPROGRAMS  CALLED: 


INPUT: 


TABX  array 
TABY  array 
TAB 2  array 


OUTPUT: 


f  |  iu*ww*»  ****-*»<*-■ 


SUBROUTINE  DISCOT 

Single  or  double  interpolation  for  continuous  or 
discontinuous  function  using  Lagrange's  formula 

CALL  DISCOT  (XA,  2A,  TABX,  TABY,  TABZ,  NC,  NY,  NZ, 
ANS) 

UNS,  DISSER,  LAGRAN 

These  subroutines  are  concerned  with  the  interpo¬ 
lation,  and  are  not  documented  separately 


x  value  (first  independent  variable)  for  interpolated 
point 

z  value  (second  independent  variable)  for  interpolated 
point 

Some  as  x  value  for  single-line  function  interpolation 

Table  of  x  values--first  independent  variable 

Table  of  y  values— dependent  variable 

Table  of  2  values — second  independent  variable 

Three  digit  control  integer  with  +  sign 

Use  +  sign  if  NX  ■  NY/NZ  *  points  in  X  array 

Use  -  sign  if  NX  -  NY 

Use  1  In  hundreds  position  for  no  extrapolation 
above  maximum  Z 

Use  0  in  hundreds  position  for  extrapolation 
above  maximum  Z 

Use  1-7  in  tens  position  for  degree  of  interpolation 
deoired  in  X  direction 

Use  1-7  in  units  position  for  degree  of  interpolation 
desired  in  Z  direction 

Number  of  points  in  y  array 

Number  of  points  in  z  array 


y  value  (dependent  variable)  interpolated  at  x,  z 

DISCOT  is  a  "standard"  routine  used  at  DTNSRDC. 
Consult  User  Services  Branch  of  the  Computation, 
Mathematics  and  Logistics  Department  for  additional 
information. 
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NAME:  FUNCTION  MINP 

PURPOSE:  Select  index  of  minimum  x  value  to  b«  used  for 

Lagrange  interpolation,  from  an  array  of  x  values 
greater  than  required 

CALLING  SEQUENCE:  I  -  MINP  (M,  N,  XA,  X) 

INPUT: 

M  m  "  number  of  points  required  for  interpolation  of 

degree  m-1 

N  n  -  total  number  of  points  in  x  array  >  m 

XA  x  value  to  be  used  for  interpolation 

X  array  Table  of  x  values,  must  be  in  ascending  order,  but 

need  not  be  equally  spaced 

OUTPUT: 

MINP  Index  of  minimum  x  value  from  the  array  to  be  used  by 

FUNCTION  YINTE  for  Lagrange  interpolation  of 
degree  m-1 

SAMPLE  PROGRAM  USING  FUNCTIONS  MINP  AND  YINTE: 

DIMENSION  X(10),  Y(10) 

N  -  10 
M  «  4 

READ  (5,  10)  (X(J),  >1,  N),  (Y(J),  J-l,  N) ,  XA 
I  -  MINP  (M,  N,  XA,  X) 

YA  -  YINTE  (XA,  X(I),  Y(I),  M) 


ALTERNATE  PROGRAM  USING  FUNCTION  YINTX: 

DIMENSION  X(10),  Y(10) 

N  -  10 
M  -  4 

READ  (5,  10)  (X(J),  J-l,  N),  (Y(J),  J-l,  N) ,  XA 
YA  -  YINTX  (XA,  X,  Y,  M,  N) 

The  result  from  either  program  is  the  same.  In  either  case,  only  the 
M  points  closest  to  XA  are  considered  in  the  interpolation  formula.  The 
first  combination  should  be  used  whenever  several  dependent  variables  are 
to  be  interpolated  at  some  value  of  the  independent  variable,  since  MINP 
need  only  be  called  once,  FUNCTION  YINTE  may  be  used  alone  whenever 
N  -  M. 
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NAME: 

PURPOSE: 

CALLING  SEQUENCE: 
INPUT: 

XA 

X  array 

Y  array 
N 

OUTPUT: 

YINTE 


'•  .1.1 


FUNCTION  YINTE 

Single  interpolation  of  degree  n-1  for  function 
represented  by  n  (x,y)  points  using  Lagrange's 
formula 

YA  -  YINTE  (XA,  X,  Y,  N) 

x  value  (independent  variable)  for  interpolated  point 

Table  of  x  values — independent  variable 

x  values  can  be  in  either  ascending  or  descending 
order  and  do  not  need  to  be  equally  spaced 

Table  of  y  values — dependent  variable 

n  *  number  of  (x,y)  values  defining  the  function 


Interpolated  y  value  (dependent  variable)  derived 
from  Lagrange  formula  of  degree  n-1 

For  example,  when  n  -  4,  cubic  interpolation  is 
performed 


Lagrange's  Interpolation  Formula 


(x-Xj)  (x~x2) 


U-Kn> 


(Xq-j^)  (x0-x2)  .  .  .  (x„-xj 


(x— Xq)  (x-x2)  . 


0*n' 


(xrV  (xrx2>  •  • 


.  (xrxn)  yl 


(x-xQ)  (x-x1)  (x-x3)  .  .  .  (x-xn) 

+  (x2-xQ)  (x2-x1)  (x2-x3)  .  .  .  (x2“xn^ 


(x-xQ)  (x-x^  (x-x2)  .  .  .  (x-x^) 

+  (xn-x0nVxl>  (xn~x2^  ’  •  “  <VVl) 
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NAME: 


FUNCTION  YINTX 


PURPOSE:  Single  interpolation  of  degree  m-1  for  function 

represented  by  n  (x,y)  points  using  Lagrange's 
formula.  If  n  >  m,  only  the  m  closest  points  are 
considered  In  the  interpolation  formula 

CALLING  SEQUENCE:  YA  -  YINTX  (XA,  X,  Y,  M,  N) 

INPUT: 

XA  x  value  (Independent  variable)  for  interpolated  point 

X  array  Table  of  x  values-- Independent  variable  x  values  must 

be  in  ascending  order,  but  need  not  be  equally  spaced 

Y  array  Table  of  y  values — dependent  variable 

M  m  *  number  of  (x,y)  values  considered  for  the 

interpolation  process  of  degree  m-1 

N  n  ■  total  number  of  (x,y)  values  >  m 

OUTPUT: 

YINTX  Interpolated  y  value  (dependent  variable)  derived 

from  Lagrange  formula  of  degree  m-1 

FUNCTION  YINTX  may  be  used  instead  of  FUNCTION  MINP 
and  FUNCTION  YINTE  together 

See  Sample  Programs  using  these  three  functions 
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APPENDIX  B 


SAMPLE  INPUT  AND  OUTPUT  FOR  PHPRIM 
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PLANtNa  HUM.  WITH  PROPELLERS  PLANING  CRAFT  WITH  2  PROPELLERS .^D  DIESEL  CNQINCS 
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SAMPLE  OUTPUT  PAGE  10 


INITIAL  DISTRIBUTION 


Copies 

14  NAVSEA 

2  PMS  300T 
1  941 

1  521 
1  5213 

1  5241 

1  03D3 
1  03D4 

1  03F2 

1  03R 

1  03R3 

1  312 
1  3124 

1  31241 


50 

NAVSEADET  Norfolk/6660 

12 

DDC 

8 

ABC-17 

CENTER  DISTRIBUTION 

Copies 

Code 

1 

11  (Ellsworth) 

1 

1103  Data  Bank 

1 

112 

1 

117 

1 

118 

1 

15 

1 

152 

1 

1524 

7 

1524  (Hubble) 

1 

1532 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY.  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


